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ABSTRACT 

We present an analysis of '^Li/^Li isotope ratios along four sight lines that probe diffuse molecular 
gas near the supernova remnant IC 443. Recent gamma-ray observations have revealed the presence 
of shock-accelerated cosmic rays interacting with the molecular cloud surrounding the remnant. Our 
results indicate that the ^Li/^Li ratio is lower in regions more strongly affected by these interactions, 
a sign of recent Li production by cosmic rays. We find that ^Li/^Li « 7 toward HD 254755, which is 
located just outside the visible edge of IC 443, while ^Li/^Li k, 3 along the line of sight to HD 43582, 
which probes the interior region of the supernova remnant. No evidence of ^Li synthesis by neutrino- 
induced spallation is found in material presumably contaminated by the ejecta of a core-collapse 
supernova. The lack of a neutrino signature in the ^Li/^Li ratios near IC 443 is consistent with recent 
models of Galactic chemical evolution, which suggest that the i^-process plays only a minor role in Li 
production. 

Subject headings: ISM: abundances — ISM; atoms — ISM: individual objects (IC 443) — ISM: su- 
pernova remnants 



1. INTRODUCTION 

The production of Li involves contributions from a 
number of different astrophysical sources, including the 
Big Bang, Galactic cosmic rays (GCRs), red giant branch 
(RGB) and asymptotic giant branch (AGB) stars, and 
Type II supernovae (SNe II). The less abundant isotope, 
^Li, is mainly a product of spallation (and a -\- a fu- 
sion) reactions induced by GCRs in interstell ar gas (e.g., 
IMeneguzzi et al. IHQTH IRamatv et al. II1997D . Evidence 
for this can be seen in the solar system abundance ra- 
tio of ^Li to ^Be, the sole stable isotope of Be which 
can be produced only thro ugh GCR spallation. The me- 
teoritic ^Li/^Be ratio (5.6; iLodders! 2003) is essentially 
the GCR value. In contrast, the ^Li/^Li ratio for solar 
system material is ~12, while standard GCR nucleosyn- 
thesis predicts a ratio of ~1.5. Since the primordial ^Li 
abundance arising from the Big Bang is approximately 
a factor of 10 less than the present-day abundance, the 
majority of Li must be produced in stars. However, the 
precise nature of the stellar source remains unclear. 

The most promising candidates for a stellar ^Li 
source include RGB and AGB stars, where ^Li is 
produced via the Cameron-F owler mechanism (e.g., 
iSackmann fc Boothrovd1ll999[) . and SNe II, where ^Li 
is synthesized by neutrino-induced spallation in the He 
and C shells of the prog enitor star during co re col- 
lapse (i.e., the I'-process; iWooslev et al. I I1990D . Di- 
rect observational evidence for the i^-process remains 
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elusive, yet could be revealed by detailed studies of 
''Li/^Li isotope ratios in interstellar gas surrounding su- 
pernova remnants (SNRs). Since virtually no ^Li is ex- 
pected to be produced by neutrino nucleos ynthe sis (e.g., 
IWooslev fc Weaver I [l995l : lYoshida et am2008l ). inter- 
stellar material contaminated by SN II ejecta should ex- 
hibit a ''Li/^Li ratio that is enhanced over the ambient in- 
terstellar value. However, a reduced ^Li/^Li ratio might 
also be anticipated near a supernova remnant because 
SNRs are thought to be the primary sources responsible 
for cosmic-ray acceleration. Interactions between cosmic 
rays, accelerated by a supernova shock, and any nearby 
interstellar gas will drive the ^Li/^Li ratio toward the 
pure GCR value. 

The well-studied SNR IC 443 provides an excel- 
lent opportunity to test which of these processes dom- 
inates the production of Li by SNe II. Located at 
a distance of 1.5 kpc in the Gem OBI associa- 
tion, IC 443 is an intermediate-age (10—30 kyr) core- 
collapse SNR known to be interacting with atomic 
and molecular gas in its vicinity. Evidence of the 
interaction in cludes observations of shocked H i fil- 
aments (e.g., iBraun fc Strom I Il986t iLee et al. I l2008fl 
and shocked molecula r clumps fe.g..lHuang et al. II1986I : 
lDickmanet"anil992t iSnell et al. 1 120051 ). as well as the 
detection of OH (1720 MHz) masers coincident with 
the shocked molecu lar material (jClaussen et al. 1 119971 : 
iHewitt et al. I 120061 ). Many recent investigations of 
IC 443 have focused on the emis sion of gamma rays 
at GeV and TeV energies (e.g . . lAcciari et ai~1 120091 : 



iTavani et alTI [20101: lAbdo et al. I |2010D . These studies 



strongly suggest a pionic origin for the gamma-ray emis- 
sion, indicating that shock-accelerated cosmic rays (in 
addition to the shocks themselves) are interacting with 
the ambient molecular cloud. 

In this Letter, we present an analysis of ^Li/^Li iso- 
tope ratios along four sight lines through IC 443. The 
ratios are extracted from high-resolution observations of 
Li I A6707 toward HD 254477, HD 254577, HD 43582, 
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Figure 1. Composite image of IC 443. The optical image from 
the Digitized Sky Survey (POSS-II/red filter) is shown superim- 
posed onto a gamma-ray intensity map in the 400 MeV to 3 GeV 
energy range obtained by the Fermi LAT (J. Hewitt 2011, private 
communication). The stars targeted for HET observations are la- 
beled. 

and HD 254755. Figure 1 shows the positions of the 
target stars in relation to the gamma-ray emission de- 
tected by the Fermi Large Area Telescope (LAT) in 
the 400 MeV to 3 GeV energy range (J. Hewitt 2011, 
private communication). Previous studies of the opti- 
cal absorption profiles toward the sel ected targets have 
revealed extensive high- velo city gas () Welsh fc Sallmen"! 
I2003t iHirschauer et al. II2009D , indicating that these stars 
probe material directly associated with the SNR. 

2. OBSERVATIONS AND DATA REDUCTION 

The target stars were observed with the 9.2 m Hobby- 
Eberly Telescope (HET) at McDonald Observatory dur- 
ing the winters of 2008/2009 and 2009/2010. AU ob- 
servations employed the High Resolution Spectrograph 
ITuII 1998) with an effective slit width of 125 //m 
{R ss 120, 000) so that the fine structure lines of Li i 
could be adequately resolved. Two spectrographic set- 
tings (centered at 4931 A and 5936 A) provided data 
on Li I A6707. The setting at shorter wavelengths was 
needed to simultaneously obtain data on K i A4044, while 
the longer wavelength setting allowed the K i A7698 line 
to be observed. Both settings also yielded information 
on Ca I A4226, CH+ A4232, and CH A4300. 

Basic i nformation concer n ing th e target stars can be 
found in IHirschauer et al.~l (|2009D . The three brighter 
targets were observed with both instrumental setups, re- 
sulting in total exposure times on Li i A6707 of 6.1 hr 
for HD 254577 and 4.7 hr for HD 43582 and HD 254755. 
These exposure times yielded signal-to-noise ratios (S/N) 
per resolution element near Li i of about 1100. The 
faintest star in our sample, HD 254477, was observed 
only with the longer wavelength setting for a total of 2.0 



hr, resulting in a S/N of 380. 

Standard procedures within IRAF were employed for 
bias correction, cosmic-ray removal, scattered light sub- 
traction, and flat fielding. Ideally, a master flat was ob- 
tained from the median of flats taken for a given night. 
In some instances, only the first flat of a sequence could 
be used due to the appearance of an emission feature 
near 6708 A in subsequent flats. One-dimensional spec- 
tra were extracted from the processed images and were 
wavelength-calibrated after identifying emission lines in 
Th-Ar comparison spectra. Before shifting the calibrated 
spectra to the local standard of rest (LSR) frame, small 
velocity corrections were applied based on the measured 
wavelengths of the atmospheric [O i] emission lines at 
6300 A and 5577 A. This was necessary to account for 
offsets in the velocity zero points caused by slight changes 
in spectrograph placement between the stellar and Th-Ar 
exposures. 

Final spectra were produced by co-adding the indi- 
vidual exposures of a given target for orders containing 
the lines of interest. The co-added spectra were normal- 
ized to unity by fitting low-order polynomials to regions 
free of interstellar or telluric features. Figure 2 presents 
the reduced Li i spectra for the four sight lines probing 
IC 443. From the widths of Th i emission lines in the 
nightly comparison spectra, we found that the actual re- 
solving power of the spectrograph during our observa- 
tions was R = 98, 000. 

3. ANALYSIS 

Detailed knowledge of the velocity component struc- 
ture in each direction is necessary if meaningful ^Li/^Li 
ratios are to be extracted from the complicated Li I line 
profiles, which exhibit fine and hyperfine structure in 
addition to isotopic splitting. For a single, optically- 
thin cloud, and in the absence of ^Li, the two fine- 
structure components of the Li i doublet will exhibit rela- 
tive strengths of 2:1 and a separation of 0.15 A. Addition 
of ^Li reduces this ratio because the stronger of the two 
^Li components is superimposed onto the weaker of the 
components from ^Li (i.e., the isotope shift is approxi- 
mately equal to the fine-structure separation). 

Further complications arise when multiple clouds are 
present along the line of sight. Thus, the first step in 
our analysis was to derive a robust solution for the line- 
of-sight component structure using species with mod- 
erately strong absorption lines. The K i A7698 and 
CH A4300 lines were of particular interest as these species 
ar e expected to coexist with Li i in cool, di f fuse clouds 
(iWeltv fc Hobbs l[200l [Knauth et al~l[200l |Pan et al. I 
I2005f) . While the weaker K i line at 4044 A would provide 
a better template for Li i A6707 (due to these lines having 
similar intrinsic strengths), absorption from K i A4044 
is detected only toward HD 254755. We also analyzed 
the Ca I A4226 and CH+ A4232 absorption profiles to 
check for consistency in velocity among the components 
in common with K I and CH. 

The column densities and component parameters for 
all of the relevant species were determined through pro- 
file synthesis using the program ISM OD developed by 
Y. Sheffer (see, e.g.. ISheffer etaTIIMlS) . ISMOD treats 
the velocities, 6-values, and column densities of the ab- 
sorption components as free parameters while minimizing 
the rms deviations in the residuals of the fit. Hyperfine 
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Figure 2. Profile synthesis fits to the fine structure lines of '^Li I and "^Li I toward HD 254755, HD 43582, HD 254577, and HD 254477. 
Along with the best-fit synthetic profile (black line), two alternative syntheses are shown in each case, one assuming a solar system ^Li/®Li 
ratio (red line) and one assuming no ®Li is present (grey line). The parameters for ^Li in these alternative syntheses are the same as in 
the best-fit cases. Residuals are given above each fit. For HD 43582 and HD 254577, the best-fit profile displayed is that based on the K I 
template. Note that the alternative syntheses are less satisfactory in general, but particularly for HD 43582. 



structure is included in our synthetic profiles for Li i 
and K i , with wavelengths and /-values adopted from 
IMorton I (|2003l ). Our fits also account for A-doubling in 
the CH line. To extract the ^Li/^Li ratio, the profile syn- 
thesis routine performs a simultaneous fit to the blended 
fine structure lines of ''Li i and ^Li i. The component 
structure is assumed to be identical for the two isotopes, 
with the isotope ratio left as a free parameter. Our re- 
sults for K 1 and CH were used to determine the number 
of velocity components to include in the synthesis for Li i. 

For HD 254755, a single component contributes more 
than 90% of the total column density of K i and CH. 
Our best one-component fit to the Li i profile (Figure 2) 
yields a ^Li/^Li ratio of 7.1 ± 2.4, where the errors are 
dominated by the observational uncertainties associated 
with the detection of ^Li. The velocity and &-value re- 
sulting from the fit (see Table 1) are in good agree- 
ment with those for the corresponding component in 
K I (wLSR = -6.4 km s'^ b = 1.2 km s^i) and CH 
(wLSR = -5.9 km s-^ b = 2.1 km s^^). 

For HD 43582 and HD 254577, two velocity compo- 
nents comprise 70—80% of the total columns of K i and 
CH. The Li i profiles in these directions are manifestly 
more blended than toward HD 254755, indicating that 
multiple components are present. We therefore kept the 
component structure fixed in the Li i syntheses, with 
component parameters determined separately from the 
two dominant components in K i and CH. The results are 
given in Table 1 and details concerning the fits are pre- 



sented in Figures 3 and 4 for HD 43582 and HD 254577, 
respectively. In both cases, the ^Li/®Li ratios derived 
from the K i and CH templates are indistinguishable 
given the uncertainties. Thus, we adopt the average 
ratios for further analysis (i.e., ^Li/^Li — 3.1 ± 1.4 for 
HD 43582 and 6.1 ± 3.2 for HD 254577). 

For HD 254477, the strongest velocity component ac- 
counts for only about 60% of the K i and CH col- 
umn density, but no other component makes up more 
than 10% of the profile in these species. Our best one- 
component fit to the Li 1 profile yields a consistent ve- 
locity and fe-value compared to the dominant compo- 
nent in K I (ulsr = —2.8 km s~^; b = 2.3 km s^^) 
and CH (wLSR = -2.5 km s'^; b = 1.0 km s^i). 
However, the significantly lower S/N makes it impossi- 
ble to confidently derive a ^Li/^Li ratio in this direc- 
tion. Nominally, we find "^Li/^Li = 4.5 ± 3.6, but with 
N{^U i) = (2.2 ± 1.8) X 10^ cm-2 (i.e., not even a 2a de- 
tection). A more appropriate 3a upper limit on A^(^Li i) 
of < 5.4 X 10^ cm~^ gives a lower limit on ^Li/°Li of 
>l-9. 

4. DISCUSSION AND CONCLUSIONS 

Our results on ^Li/^Li isotope ratios in the vicinity 
of IC 443 shed new light on Li production associated 
with SNe II. Before these observations, it was not clear 
whether one should expect an enhanced abundance of 
^Li near SNe II, resulting from neutrino spallation during 
the core-collapse phase, or a higher relative abundance 
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Figure 3. Synthesis fits to the Li I profile toward HD 43582, 
based on the K I template (upper panel) and the CH template 
(lower panel). The dot-dashed lines show the contributions from 
individual cloud components to the overall profile, shown as a solid 
line. 

of ^Li, due to spallation and fusion reactions initiated 
by shock-accelerated cosmic rays. Our analysis of Li i 
absorption along four sight lines through IC 443 suggests 
that cosmic-ray interactions dominate the production of 
Li at least for this particular SNR. 

All of the ^Li/^Li ratios we obtain near IC 443 are 
lower than the solar sys tem (meteoritic ) valti e of 12.2 
(jLodders I 2003), which IKnauth et all (|2003D suggest 
is representative of gas in the solar neighborhood. 
While our results for HD 254577 and HD 254755 
are actually consistent with the solar system ratio at 
the la level, the ratio we find toward HD 43582 is 
lower by more than 6(j. There is evidence, however, 
that the Li isotope ratio in the interstellar medium 
(ISM) in general is lower than the solar system value. 
The weighted mean ISM rati o from measur e ments 
available in the literature (iLemione et al. I 119931: 
Mever et al. 1 119931: iLemoine et al 1 119951: IKnauth et al. I 



2000"; 'Howarth et al. I 120021: IKnauth et al. I 120031: 
Kawanomoto et al. I I2009I) is 7.3 ± 0.6 (if one considers 
only the strongest absorption component in each case). 
This value is very similar to the ratio we find toward 
HD 254755, which lies beyond the outer edge of the 
Ha emission contours in the northeast of IC 443 (see 
Figure 1). The gas in this direction may thus represent 
material relatively unaffected by the supernova shock, 
implying that a ^Li/^Li ratio of ^1 characterized the 
ambient molecular cloud before the occurrence of the 
supernova. 

In contrast, the ^Li/^Li ratio toward HD 43582 is 
40—80% lower than the mean value observed in the local 
ISM. The stars HD 254755 and HD 43582 are separated 




0.990 



0.985 



(CH Fit) 



6707.2 6707.4 6707.6 6707.8 6708.0 6708.2 
Wavelength (A) 

Figure 4. Same as Figure 3 except for the Li I profile toward 
HD 254577. 



by only 7'. 4 on the sky (or 3.2 pc at the distance of 
IC 443), but the line of sight to HD 43582 clearly pene- 
trates the interior region of the SNR. I n fact, HD 43582 is 
the only star of those investigated by iHirschauer et al.l 
(2009) that exhibits Ca ii absorption at high positive 
velocity, strongly suggesting that it lies behind the re- 
ceding edge of the expanding shell of the remnant. If 
the molecular gas in this direction has been subjected 
to a more intense cosmic-ray irradiation than that to- 
ward HD 254755, then an enhancement in the relative 
abundance of ®Li would be a natural consequence. In- 
deed, the ^Li/^Li ratio toward HD 43582 approaches 
the value predicted by models of GCR nucleosynthe- 
sis (^1.5: i Mencgu zzi et al. 111973: iRamatv et al. 1 119971: 
ILemoine et al. 199^). The sight lines to HD 254577 and 
HD 254477 presumably also pass through the interior of 
the SNR. However, given the uncertainties, the ^Li/^Li 
ratios in these directions are consistent with both the 
mean ISM value (as observed toward HD 254755) and 
the lower value we find toward HD 43582. 

The discovery of a low ^Li/^Li ratio in IC 443 cor- 
roborates the conclusions drawn fr om gamma-ray emis- 
sion studies of th e region (e.g . . lAcciari et al. I 120091 : 



iTavani et ahl 120101 : lAbdo et all |2010D , which~strongly 
suggest that the gamma radiation results from the de- 
cay of neutral pious produced through cosmic-ray in- 
teractions with molecular gas. Complementary results 
on the co smic-ray ionization rate were recently reported 
bv I ndriolo et al. I (^OljD). These authors find large H^ 
column densities near IC 443 and deduce an H2 ioniza- 
tion rate of ^2 ~ 2 x 10"""^^ s""'^, or about five times the 
rate typically found for diffuse molecular clouds (e.g., 
Ilndriolo et aUl2007l) . 

Both the HJ abundances discussed bv Ilndriolo et al. I 
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(|201ClD and the ^Li/^Li ratios presented here indicate 
regions of enhanced cosmic-ray activity within IC 443. 
However, the two observables present potentially inter- 
esting differences in term s of which sight line s exhibit 
the largest enhancements. Ilndriolo et al. I (J201C0 infer an 
ionization rate of C2 = 2.6iy^xlCr^s^ (the highest in 
their survey) toward HD 254577, yet report only an up- 
per limit on C2 toward HD 43582 since H3 is not detected. 
(The upper limit is < 1.4 x 10"^^ s~^ if uncertainties in 
the H2 column density are included.) These results ap- 
pear to be in conflict with our determinations of ^Li/^Li 
for these directions, since we find a lower ^Li/^Li ratio 
toward HD 43582 than toward HD 254577. However, the 
uncertainties both in our determinations and in those of 
Indriolo et al., which rely on a number of assumptions 
concerning the physical conditions in the clouds contain- 
ing H;^, are large. Moreover, since the abundance of H;^ 
traces the current cosmic-ray ionization rate, while the Li 
isotope ratio probes the integrated cosmic-ray flux, some 
differences in these two observables might be expected 
(B. Fields 2011, private communication). 

A decrease in the ^Li/^Li ratio due to enhanced cosmic- 
ray activity (resulting in newly-synthesized ^Li and ''Li) 
should be accompanied by an increase in the elemental Li 
abundance. If the clouds toward HD 43582 initially pos- 
sessed a "^Li/^Li ratio of ~7, then the presently-observed 
ratio of 3.1 would imply a factor-of-two increase in the el- 
emental Li abundance. Some evidence for this effect can 
be seen in the A^(Li i)/A^(K i) ratios found toward the 
stars in our sample: (2.7 ± 0.2) x 10"^ ^j. jjD 254755, 
(5.9 ± 0.6) X 10-3 for HD 254577, (6.4 ± 0.9) x lO'^ 
for HD 43582, and (6.8 ± 1.4) x IQ-^ for HD 254477. 
The three sight lines that pass through the interior of 
the remnant seem to exhibit the expected factor-of-two 
enhancement in the abundance of Li (relative to K) if 
the initial abundances were similar to that seen toward 
HD 254755. 

Still, none of the A^(Li i)/A^(K i) ratios we find near 
IC 443, are entirely unusual compared to value s typical of 
diffus e clouds (e.g., Welty & Hobbs 2001; Kn auth et al. I 
120031 . Without directly determining the elemental Li/K 
ratios, which would require making uncertain corrections 
for ionization and depletion onto grains, it is difficult to 
conclude whether an enhancement in Li is indeed ob- 
served. We can, however, state that a localized enhance- 
ment would be energetically feasible. If we assume that 
riH ^ 200 cm~'^ for the g as toward HD 43582 (as adopted 
by Ilndriolo et al. 1 120101) , then a factor-of-two increase 
in the Li abundance would require that ~3 x 10^° erg 
are impart ed to the accelerated cosmic rays by the SN 
shock fsee iRamatv et al. |[2000l) . This is only -20% of 
the typical mechanical en ergy available in SN II ejecta 
(|Wooslev fc Weaver |[T995h 

Our observations suggest that the abundance of ^Li 
relative to ^Li has been enhanced near IC 443 by inter- 
actions between shock-accelerated cosmic rays and the 
ambient molecular cloud. At the same time, we find 
no evidence of ^Li synthesis by neutrino-induced spal- 
lation in material that presumably has been contami- 
nated by the ejecta of a core-collapse supernova. Given 
the age of the remnant, the hot ejecta should have had 
ample time to cool, as it interacts with its dense sur- 
roundings, to the point at which absorption lines from 
neutral atoms would become visible. For comparison. 



iWallerstein et al.1 (|1995f ) estimated a cooling time of 800 
yr for the Vela remnant, which is about the same age 
as IC 443 but where the ambient ISM is of much lower 
density. The lack of a neutrino signature is consistent 
with recent mod els of Galactic chemical evolution (e.g., 
IPrantzosI I2012D . which suggest a very minor role for 
the i^-process in producing ^Li. Future measurements 
of ^Li/^Li in gas surrounding other SNRs will help to 
establish unequivocally the role that neutrino spallation 
plays in Li production. 
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Table 1 

Li I Profile Synthesis Results 






Star 


Template 


t'LSR 
(km s-i) 


b 

(km s-l) 


nCu i) 

(10^ cm-2) 


^Li/6Li 


Af(Li i) 
(IQS cm-2) 


HD 254477 




-2.2 


1.7 


10.2 ±1.8 


> 1-9 


12.4 ±2.5=^ 


HD 254577 


K I 
CH 


(-6.9) 
(-2.9) 
(-6.7) 
(-2.5) 


(2.3) 
(2.4) 
(1.8) 
(1.7) 


5.8 ±0.7 
4.2 ±0.7 
6.1 ±0.6 

3.9 ±0.6 


5.7±3.1 
6.4 ±3.4 


11.8 ±1.3 
11.6±1.1 


HD 43582 


K I 
CH 


(-7.4) 
(-3.9) 
(-6.9) 
(-2.9) 


(1.4) 
(0.6) 
(2.0) 
(1.8) 


4.3 ±0.6 
1.8 ±0.5 
5.1 ±0.7 
1.3 ±0.6 


2.8 ±1.0 
3.3 ±1.7 


8.3 ±1.1 
8.3 ±1.3 


HD 254755 




-6.3 


2.2 


16.6 ±0.8 


7.1 ±2.4 


18.9 ±1.1 



Note. — Velocities and fe-values shown in parentheses were held fixed during profile 
synthesis. The last column gives the total Li I column density. The corresponding K I 
column densities for components also seen in Li I are 18.3±0.1, 19.7±0.1, 12.9±0.1, and 
69.4±1.7 (in units of 10" cm-2) for HD 254477, HD 254577, HD 43582, and HD 254755, 
respectively. 

* Assumes '^Li/^Li ~ 4.5. 



